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ABSTRACT 

Context. Solar-like oscillations are stochastically excited by turbulent convection at the surface layers of the stars. 
Aims. We study the role of the surface metal abundance on the efficiency of the stochastic driving in the case of the 
CoRoT target HD 49933. 

Methods. We compute two 3D hydrodynamical simulations representative - in effective temperature and gravity - of 
the surface layers of the CoRoT target HD 49933, a star that is rather metal poor and significantly hotter than the 
Sun. One 3D simulation has a solar metal abundance, and the other has a surface iron-to- hydrogen, [Fe/H], abundance 
ten times smaller. For each 3D simulation we match an associated global ID model, and we compute the associated 
acoustic modes using a theoretical model of stochastic excitation validated in the case of the Sun and a Cen A. 
Results. The rate at which energy is supplied per unit time into the acoustic modes associated with the 3D simulation 
with [Fe/H]=-1 is found to be about three times smaller than those associated with the 3D simulation with [Fe/H]=0. 
As shown here, these differences are related to the fact that low metallicity implies surface layers with a higher mean 
density. In turn, a higher mean density favors smaller convective velocities and hence less efficient driving of the acoustic 
modes. 

Conclusions. Our result shows the importance of taking the surface metal abundance into account in the modeling of 
the mode driving by turbulent convection. A comparison with observational data is presented in a companion paper 
using seismic data obtained for the CoRoT target HD 49933. 
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1. Introduction 

Using the measured Unewidths and the ampUtudes of the 
solar acoustic modes, it has been possible to infer the rate at 
which energy is supplied per unit time into the solar acous- 
tic modes. Using these constraints, different models of mode 
excitation by turbulent convection have been extensively 
tested in the case of the Sun (see e.g. recent reviews by 
ISamadi et all (|2008bD and lHoudekl (|2006f )). Among the dif- 
ferent approa ches, we can distinguish pure theoretical a^ : 



proaches (e.g. lSamadi fc Goupi|2001l:IChaplin et al.| 200 
semi-analytical approaches (e.g. lSamadi et al.ll2003bl j) and 



pure nume rical a pproaches (e.g. iNordlund fc Stein T2OOII ; 
IStein et al.ll2004t (jacoutot et"alll2008f) . The advantage of 
a theoretical approach is that it easily allows massive com- 
putation of the mode excitation rates for a wide variety 
of stars with different fundamental parameters (e.g. effec- 
tive temperature, gravity) and different surface metal abun- 
dance. However, pure theoretical approaches are based on 
crude or simplified descriptions of turbulent convection. 
On the other hand, a semi-analytical approach is gener- 
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ally more realistic since the quantities related to turbulent 
convection are obtained from 3D hydrodynamical simula- 
tion. 3D hydrodynamical simulations are at this point in 
time too time consuming, so that a fine grid of 3D models 
with a sufRcient resolution in effective temperature (Toff), 
gravity (logg) and surface metal abundance {Z) is not yet 
available. In the present paper, we study and provide a pro- 
cedure to interpolate for any value of Z the mode excitation 
rates T between two 3D simulations with different Z but 
the same Toff and logg. With such interpolation procedure 
it is no longer required to have at our disposal a fine grid 
in Z of 3D simulations. 

The semi-analytic a l mod e that we consider here is based 
on ISamadi fc Goupill (|200lD's theoretical model wi th the 
improvements proposed by iBelkacem et al.l ([2006^ . This 
semi-analytic al model satisfactorily reproduces the sola r 
seismic data (jSamadi et al.ll2003al: iBelkacem et al.ll2006b[ ). 
Recently, the seismic constraints obtained for a Cen A 
(HD 128620) have provided an additional validation of 
the basic physical as sumptions of this theoretical model 
([Samadi et al.l l2008af ). The star a Cen A has a sur- 
face gravity (log g = 4.305) lower than that of the Sun 
(log (70 — 4.438), but its effective temperature (Toff = 
5810 K) does not significantly differ from that of the Sun 
(Tcff,0 = 5780 K). The higher Toff, the more vigorous the 
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convective velocity at the surface and the stronger the driv - 
ing by turbulent convection (see e.g. iHoudek et al.|[l999f ). 
For main sequence stars with a mass M < 1.6 Mq, an in- 
crease of the convective velocity is expected t o be associated 
with a larger turbulent Mach number, Mt (|Houdek et al.l 
Il999t ). However, the theoretical models of stochastic exci- 
tation are strictly valid in a medium where Mt is - as in the 
Sun and a Cen A - rather small. Hence, the higher Mt, the 
more questionable the different approxi mations and the as- 
sump tions involved in the theory (see e.g. lSamadi fc Goupill 
l200lh . It is therefore important to test the theory with an- 
other star characterized by a Toff significantly higher than 
in the Sun. 

Furthermore, the star a Cen A has an iron-to-hydrogen 
abunda nce slightly larger than the Sun, na mely [Fe/H] = 
0.2 (see lNeuforge-Verheecke fc MagainI 19971). However, the 
modeling performed bv lSamadi et all (^2008a^ for a Cen A 
assumes a solar iron abundance ([Fe/H]=0). According to 
iHoudek et al.l (|1999[ ) , the mode amplitudes a re expected to 
chang e with the metal abundance. However. IHoudek et al.l 
(|1999f )'s result was obtained on the basis of a mixing- length 
approach involving several free parameters and by using a 
theoretical model of stochastic excitation in which a free 
multiplicative factor is introduced in order to reproduce the 
maximum of the sola r mod e excitati on rat es. Therefore, it is 
important to extend [Houdek et al] (|l999l )'s study by using 
a more realistic modeling based on 3D hydrodynamical sim- 
ulation of the surface layers of stars and a theoretical model 
of mode driving that reproduces - without the introduction 
of free parameters - the available seismic constraints. 

To this end, the star HD 49933 is an interesting case 
for three reasons: Fir st, this star has Toff = 6780 ± 130 
K (iBruntt et al.ll2008f ). logg ~ 4.25 ± 0.13 (iBruntt et all 
20081) and "[fVHI - -0.37 dex (|Solano et al.l l2005t 



Gillon fc MagainI [200l . The properties of its surface lay- 



ers are thus significantly different from those of the Sun 
and a Cen A. Second, HD 49933 was observed in Doppler 
velocity with the HARPS spe ctrograph. A seismi c analysis 
of these data performed by iMosser et al.l (|2005[ ) has pro- 
vided the maximum of the mode surface velocity (l^nax)- 
Third, the star was more recently observed continuously 
in intensity by CoRoT during 62 days. Apart from ob- 
servations for the Sun, this is the longest seismic observa- 
tion ever peformed both from the ground and from space. 
This long term and continuous observation provides a very 
high frequency resolution ('^ 0.19 //Hz) . The seismic analy- 
sis of these observations u ndertaken bvlAppourchaux et al.l 
(|2008t ) or more recently bv lBenomar et all (|2009f )" have pro- 
vided the direct measurements of the mode amplitudes 
and the mode linewidths with an accuracy not previously 
achieved for a star other than the Sun. 

We consider two 3D hydrodynamical simulations rep- 
resentative - in effective temperature and gravity - of 
the surface layers of HD 49933. One 3D simulation has 
[Fe/H]=0, while the second has [Fe/H]= -1. For each 3D 
simulation, we match an associated global ID model and 
compute the associated acoustic modes and mode excita- 
tion rates, T. This permits us to quantify the variation of 
y induced by a change of the surface metal abundance Z. 
From these two sets of calculation, we then deduce T for 
HD 49933 by taking into account the observed iron abun- 
dance of the star (i.e . [Fe/H]=:-0.37). In a companion paper 
(| Samadi et al.|[2009l . hereafter Paper II), we will use these 
theoretical calculations of T and the mode linewidths ob- 



tained from the seismic analysis of HD 49933 performed 
with the CoRoT data to derive the expected mode ampli- 
tudes in HD 49933. These computed mode amplitudes will 
then be compared with the observed ones. This compari- 
son will then constitute a test of the stochastic excitation 
model with a star significantly different from the Sun and 
a Cen A. It will also constitute a test of the procedure pro- 
posed here for deriving T for any value of Z between two 
3D simulations with different Z. 

The present paper is organised as follows: we first de- 
scribe in Sect. [H the method to compute the theoretical 
mode excitation rates associated with the two 3D hydro- 
dynamical simulations. Next, the effects on J* of a different 
surface metal abundance are presented in Sect.[31 Then, by 
taking into account the actual iron abundance of HD 49933, 
we derive theoretical values of 7 expected for HD 49933. 
Finally, Sect. [5] is dedicated to our conclusions. 

2. Calculation of mode excitation rates 

2.1. Model of stochastic excitation 

The energy injecte d into a mode per unit time J* is given 
by the relation fsee lSamadi fc GouDilll200ll : lBelkacem et al.l 
l2006bD : 

y=±{c], + cl), (1) 

where C|. and Cg are the turbulent Reynolds stress and 
entropy contributions, respectively, and 



I'M 

1= / dm|^, 

"'0 



(2) 



is the mode inertia, is the adiabatic radial mode dis- 
placement and M is the mass of the star. The expressions 
for and C| are given for a radial mode with frequency 
Wosc by 



C-R = I am.— — Jr^R(r,^osc) , 
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2 167r^ f {as su)^ 
C-s = / dm^yg gr Ss{r,uJosc) 

where we have defined the "source functions" : 

SR{r,uJosc) = J -j^E (k) 

X y"dwXfc(tJ + cJosc) Xfc(w) 

w' J fc^ 

X / dwxfc(^ + ^osc)Xfc(t^) 



(3) 
(4) 



(5) 



(6) 



where P is the gas pressure, p the density, s the entropy, 
p the equilibrium density profile, as = {dP/ds)p, fr = 
(d^r/dr)^ and g^. are two functions that involve the first 
and second derivatives of respectively, k is the wavenum- 
ber, E{k) is the turbulent kinetic energy spectrum, Es{k) 
is the spectrum associated with the entropy fluctuations 
(s), s is the rms of s, Xfc is the time-correlation function 
associated with the velocity, w is a characteristic velocity 
defined in a way that 3u^ ~ {u^), {■) refers to horizontal 
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and time average, u is the turbulent veloc ity field, and fi- 
nally JC^, = is the Kurtosis fsee iBelkacem et "all 

[2006a . b.. for details). Furthermore, we have introduced for 
convenience the characteristic frequency too and the char- 
acteristic wavenumber kg: 

uq = kou (7) 
fco ^ ^ (8) 

where A is a characteris t ic size derived from E{k) as ex- 
plained in ISamadi et al.l (|2003b| ) . Note that the introduc- 
tion of the term fcj] in the RHS of Eq. ([5]) and the 
term /Cg luq s^^ in the RHS of Eq. ^ ensure dimension- 
less source functions. 

The kinetic spectr um E(k) is der i ved fro m the 3D sim- 
ulation as de tailled in lSamadi et al.l (j2003b[ ). As shown by 
ISamadi et al.i (.2003bi) , the fc-dependence of Es{k) is similar 
to th at of the Ejk) Accord ingly, we assume Eg (x E. 

In ISamadi et al.l (|2008al ) , two different analytical func- 
tions for Xki^) have been considered, namely a Lorentzian 
function and a Gaussian one. In the present study we will 
in addition d erive Xki^) directly fr om the 3D simulations 
as detailled in lSamadi et al.l ([2003^. Once Xfc(w) is derived 
from the 3D simulation, it is implemented in Eq. ([5]) and 
Eq. dSD. 

We compute th e mode excitation as detailled in 
ISamadi et al.l (|2008al ): all required quantities - except ^r, 
I and Wosc ~ are obtained directly from two 3D hydro- 
dynamical simulations representative of the outer layers of 
HD 49933, whose characteristics are described in Sect. 12.21 
below. 

The quantities related to the modes (wosc, I 
and ^r) are calculated using the adiabatic pulsation 
code ADIPLS (IChristensen-Dalsgaard fc Berthomieulll99lD 
from ID global models. The outer layers of these ID models 
are derived from the 3D simulation as described in Sect. 12.31 

2.2. The 3D simulations 

We computed two 3D radiation-hydro dynamical model at- 
mospheres with the cod e CO^BOLD (jFrevtag et al.ll2002t 
IWedemever et all |2004| ) . One 3D simulation had a solar 
iron-to-hydrogen [Fe/H]=0.0 while the other had [Fc/H]=- 
1.0. The 3D model with [Fe/H]=0 (resp. [Fe/H]=-1) will 
be hereafter referred to as model SO (resp. SI). The as- 
sumed chemical composition is similar (in particular for 
the CNO elemen ts) to that of t he solar chemical composi- 
tion proposed bv lAsplund et al] (p005i ). The abundances of 
the a-elements in model SI were assumed to be enhanced 
by 0.4 dex. For SO we obtain Z/X = 0.01830 and Y=0.249, 
and for SI Z/X = 0.0036765 and y=0.252. Both 3D sim- 
ulations have exactly the same gravity (log 5 =4.25) and 
are very close in effective temperature (Tofi). Both models 
employ a spatial mesh with 140 x 140 x 150 grid points, and 
a physical extent of the computational box of 16.4 x 16.4 x 
24.2 Mm^. The equation of state takes into account the 
ionisation of hydrogen and helium as well as the formation 
of H2 molecules according to the Saha-Boltzmann statis- 
tics. The wavelength dependence of the r adiative transfer 
is tr eated by the opacity binnin g method (jNordlundlll982t 
iLudw ie 199a lVogler etani200l using five wavelength bins 
for model SO and six for model SI. Detailed wavelength- 
dependent opacities were obtained from the MARCS model 



atmosphere package (jGustafsson et al.ll2008D . Table [T] sum- 
marizes the characteristics of the 3D models. The effective 
temperature and surface gravity correspond to the param- 
eters of HD 49933 within the observational uncertainties, 
while the two metallicities bracket the observed value. 

For each 3D simulation, two time series were built. One 
has a long duration (38h and 20h for SO and SI, respec- 
tively) and a low sampling frequency (10 mn). This time 
series is used to compute time averaged quantities (p, E{k), 
etc.). The second time series is shorter (8.8h and 6.8h for 
SO and SI, respectively), but has a high sampling frequency 
(1 mn). Such high sampling frequency is required for the 
calculation of Xfe(aj). Indeed, the modes we arc looking at 
lie between ly « 1.25 mHz and v « 2.4 mHz. 

Label [Fe/H] Y Z Z/X ng [K] 

50 0.249 13.5 10"^ 0.018305 6725 ± 17 

51 -1 0.252 2.74 10"^ 0.003676 6730 ± 12 

Table 1. Characteristics of the 3D simulations. 



The two 3D simulations extend up to T = 100 000 K. 
However, for T > 30 000 K, the 3D simulations are not 
completely realistic. First of all, the MARCS-based opaci- 
ties are provided only up to a temperature of 30 000 K; for 
higher temperatures the value at 30 000 K is assumed. Note 
that we refer to the opacity per unit mass here. For the ra- 
diative transfer the opacity per unit volume is the relevant 
quantity, i.e. the product of opacity per mass unit and den- 
sity. Since in the simulation the opacity is still multiplied at 
each position with the correct local density, the actual error 
we make when extrapolating the opacity is acceptable. 

Another limitation of the simulations is the restricted 
size of the computational box which does not allow for a 
full development of the largest flow structures, again in the 
layers above T ~ 30 000 K. Two hints make us believe that 
the size of the computational domain is not fully sufficient: 

i) in the deepest layers of the simulations there is a ten- 
dency that structures align with the computational grid; 

ii) the spatial spectral power P of scalar fields in a hor- 
izontal layer does not tend towards the expected asymp- 
totic behaviour P x k for low spatial wavenumber k. We 
noticed this shortcoming only after the completion of the 
simulation runs. To mitigate its effect in our analysis, we 
will later by default integrate the mode excitation rates 
up to T = 30 000 K. However, for comparison purposes, 
some computations have been extended down to the bot- 
tom of the 3D simulations. For SO, the layers located below 
T ~ 30 000 K contribute only by < 10 % to the excitation of 
the modes lying in the frequency range where modes have 
the most chance to be detected {1/ ~ 1.2 — 2.5 mHz). For SI, 
the contribution of the deep layers is even smaller 5 %). 

Finally, one may wonder how the treatment of the small- 
scales or the limited spatial resolution of the simulation can 
influence our calculations. Dissipativc processes are han- 
dled in CO^BOLD on the one hand side implicitely by the 
numerical scheme (Roe- type approximate Riemann solver), 
and on the other hand e xplicitely by a sub-g rid model ac- 
cording to the classical ISmagorinskvl (jl963t ) formulation. 
iJacoutot et al.l (j2008f ) found that computed mode excita- 
tion r ates signiflcan t ly de pend on the adopted sub-grid 
model. ISamadi et al.l (|2007t ) have found that solar mode ex- 
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citatio n rates computed in the manner o f iNordlund fc SteinI 
()2001f ). i.e., using data directly from the 3D simulation, de- 
crease as the spatial resolution of the solar 3D simulation 
decreases. As a conclusion the spatial resolution or the sub- 
grid model can in fluence computed mo de excitation rates 
(see a discussion in lSamadi et al.|[2b08af ) . How ever, concern- 
ing th e spatial resolution and according to I Samadi et al.l 
(|2007t )'s results, the present spatial resolution (1/140 of 
the horizontal size of the box and about 1/150 of the ver- 
tical extent of the simulation box) is high enough to ob- 
tain accurate computed energy rates. The increased spa- 
tial resolution of our models in comparison to the work of 
iJacoutot et al.l (pOOSil reduces the impact of the unresolved 
scales. 



2.3. The ID global models 
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For each 3D model we compute an associat ed ID global 
model. The models are built in the manner o f lTrampedachI 
()1997f) as detained in ISamadi et"aI1 l)2008af) in such way 
that their outer layers are replaced by the averaged 
3D simulations described in Sect. 12.21 The interior of 
the models are obtained with the CESAM code assum- 
ing st andard physics : Con vection is described accord- 
ing to iBohm-Vitensd (jl958[ )'s local mixing-length theory 
of convection (MLT), and turbulent pressure is ignored. 
Microscopic diffusion is not included. The OPAL equation 
of state is assumed. The c hemical mixture of the heavy ele- 
ments is_simil^£_tothat_of (|2005[ )'s mixture. 
As in ISamadi et al.l (|2008a[ ). we will refer to these models 
as "patched" models hereafter. 

The two models have the effective temperature and the 
gravity of the 3D simulations. One model is matched with 
SO and has [Fe/H]=0, while the second is matched with SI 
and has [Fe/H]=-1. The ID models have the same chemical 
mixture as their associated 3D simulations. The parameters 
of the ID patched models are given in Tabled The stratifi- 
cation in density and temperature of the patched ID models 
are shown in Fig. ([1]). At any given temperature the den- 
sity is larger in SI as a consequence of its lower metal abun- 
dance. Indeed, the lower the metal abundance, the lower the 
opacity ; then, at a given optical depth (r), the density is 
larger in SI compared to SO. The photosphere corresponds 
to the optical depth r = 2/3. Since the two 3D simulations 
have approximatively the same effective temperature, the 
density in SI is larger at optical depth r = 2/3. Since the 
density in SI increases with depth even more rapidly than 
in SO; the density in SI remains larger for r > 2/3 than in 
SO. 



3. Effects of the metal abundance on excitation 
rates 

The mode excitation rates (T) are computed for the two 
3D simulations according to Eqs. H])-®. The integration 
is performed from the top of the simulated domains down to 
T = 30 000 K (see Scct.[22|). In the following, (rcsp. To) 
corresponds to the mode excitation rates associated with 
the 3D model with [Fc/H]=-1 (resp. [Fe/H]=0) 



Fig. 1. Mean density p as a function of temperature, T. 
The solid line corresponds to the 3D model with the metal 
abundance (SO) and the dashed line to metal poor 3D model 
(SI). The filled dots show the location where the ID models 
have been matched to the associated 3D simulation. 

3.1. Results 

Figure[2]shows the effect of the assumed metal abundance of 
the stellar model on the mode excitation rates. Ti is found 
to be three times smaller than To, i-e. p modes associated 
with the metal poor 3D model (SI) receive approximatively 
three times less energy per unit time than those associated 
with the 3D model with the solar metal abundance (SO). 

For both 3D models, the dominant part of the driving is 
ensured by the Reynolds stresses. The entropy fiuctuations 
contribute by only 30 % of the total power for both SO 
and SI. By comparison, in the case of the Sun and a Cen A 
it contributes by only ~ 15 %. Furthermore, we find that 
the contribution of the entropy source term is - as for the 
Reynolds stress term - about three times smaller in SI than 
in SO. We conclude that the effect of the metal abundance 
on the excitation rates is almost the same for the Reynolds 
stress contribution and the entropy source term. 

3.2. Interpretation 

From Eqs. ([l]), ([2]), dS]), d?]) and (g]) we show that at a given 
layer the power supplied to the modes - per unit mass - by 
the Reynolds stress is proportional to i^kin A'* Sr/IVI, where 
-Fkin is the flux of the kinetic energy, which is proportional 
to pu^, A is a characteristic length (see Sect. 12. ip and M is 
the mode mass defined as: 



^=4 



(9) 



where is the mode displacement evaluated at the layer 
in the atmosphere where the mode is measured. 

The power supplied to the modes - per unit mass - 
by the entropy source term is proportional to p A^ 31^ Ss 
where Wosc is the mode frequency, 3? oc i^conv/^kin, where 
^conv oc w s is the convective flux, and flnally s i s the 
rms of the entropy fluctuations fsee lSamadi et al1l2006D . We 
recall that the higher 51, the higher the relative contribution 
of the entropy source to the excitation. We study below the 
role of M, Fkin, A, Sr, Ss and 3?: 
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[Fe/H] Y Z Teff [K] R/Rq M/Mq a 

0249 13.5 10"^ 6726 073 OOS 1677 
-1 0.252 2.74 10"^ 6732 1.261 1.033 1.905 

Table 2. Characteristics of the ID "patched" models, a is the mixing-length parameter. 




0.5 1.0 1.5 2.0 2.5 

u [mHz] 



Fig. 2. Mode excitation rates T as a function of the mode 
frequency, v. The solid line corresponds to the 3D model 
with the canonical metal abundance (SO) and the dashed 
line to the metal poor 3D model (SI). The dot-dashed line 
corresponds to the mode excitation rates derived for the 
specific case of HD 49933 as explained in Appendix A. 

Mode mass (JA): The frequency domain, where modes 
are strongly excited, ranges between v w 1.2 mHz and v ~ 
2.5 mHz. In this frequency domain, the mode masses M 
associated with SO are quite similar to those associated with 
SI (not shown). Consequently the differences between Ti 
and To do not arise from the (small) differences in M. 

Kinetic energy flux ('i^kinj' The larger i^kin, the larger 
the driving by the Reynolds stress. However, we find that 
the two 3D models have very similar i^kin- This is not sur- 
prising since the two 3D models have very similar effective 
temperatures. This means that the differences between Ti 
and CpQ do not arise from the (small) differences in Fkin- 

Character is tic le ngth (A): In the manner of 
ISamadi et all (I2003bl 1 we derive from the kinetic energy 
spectra E{k) of the two 3D simulations the characteristic 
length A (A = 27r/fco: see Eq. ([8])) for each layer of the 
simulated domain. We find that the differences in A 
between the two 3D simulations is small and does not 
play a significant role in the differences in CP. This can 
be understood by the fact th at SO and SI h ave th e same 
gravity. Indeed, as shown by ISamadi et al.l ((20083) - at 
a fixed effective temperature - A scales as the inverse of 
g. We conclude that the differences between CPi and To do 
not originate from the (small) differences in A. 

Source functions (Sr and Ss): The dimensionless source 
functions Sr and Ss are defined in Eqs. ([5]) and ^ re- 
spectively. Both source functions involve the eddy time- 
correlation function Xki'-^)- We define tok as the frequency 



width of Xki^)- As shown by ISamadi et al.l ([lOOSj) and 
as verified in the present case, Wfc can be eva luated as th e 
product k Uk where Uk is given by the relation (|Steinlll967t ): 

ul = dkE{k) (10) 
where E{k) is normalised as: 

r+oo 

/ dkE{k) = = 1^2 (11) 

Jo 

According to Eqs. (jlOp and (dJ), Wfc is directly proportional 
to u. At a fixed k/ko, we then have Wfc oc m fco = <^o- 

As seen above, ujq controls ujk, the frequency width of 
Xk- Then, at fixed Wosc, we can easly see from Eqs. and 
([6]) that the smaller wq, the smaller S'fl(wosc) and Ss{uJosc)- 
Since lvq = ukg = 2ttu/A and since both 3D simulations 
have approximately the same A, smaller u results directly 
in smaller loq and hence in smaller source functions. 

We have plotted in Fig. [3] the characteristic velocity u. 
This quantity is found to be up to 15 % smaller for SI com- 
pared with SO. In other words, the metal poor 3D model is 
characterized by lower convective velocities. Consequently, 
the source functions are smaller for SI compared to SO. 
Although the convective velocities differ between SO and 
SI by only 15 %, the excitation rates differ by a factor ~ 3. 
The reason for this is that he source functions, which are 
non-linear functions of m, decrease very rapidly with u. This 
is the consequence of the behavior of the eddy-time correla- 
tion Xk- Indeed, this function varies with the ratio LOosc/^k 
approximately as a Lorentzian function. This is why Xk 
varies rapidly with u (we recall that tOk oc ufco). 

In conclusion, the differences between !Pi and Tq are 
mainly due to differences in the characteristic velocity u. 
In turn, the low convective velocity in SI is a consequence 
of the larger density compared to SO. Indeed, as shown 
in Fig. [1] the density is systematically higher in SI. At 
the layer where the modes are the most excited (i.e. at 
T ~ lOOOOiC), the density is ~ 50 % higher. Since the two 
3D models have a similar kinetic energy flux (see above), 
it follows that a larger density for SI then implies lower 
convective velocities. 

Relative contribution of the entropy source term (31): 
The convective flux i^conv in SI is almost identical to that of 
SO. This is due to the fact that the two 3D simulations have 
almost the same effective temperature. Furthermore, as 
pointed out above, the differences in Fkin between SI and SO 
are small. As a consequence, the ratio Ji oc i^conv/^kin does 
not differ between the two 3D simulations. Accordingly, as 
for the Reynolds contribution, the variation of the excita- 
tion rates with the metal abundance is only due to the 
source term Ss- The latter varies with loq in the same 
manner as Sr, which is turn the reason for As a conse- 
quence, the contribution of the entropy fluctuations to show 
the same trend with the metal abundance as the Reynolds 
stress term. 
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Fig. 3. Characteristic velocity u defined in Eq. ([7]) as a 
function of temperature, T. The solid and dashed lines have 
the same meaning as in Fig. [2l The dot-dashed line corre- 
sponds to the solid line multiplied by 71, where 71 (T) = 

1/3 

{po/pi) a-nd po (resp. pi) is the mean density stratifica- 
tion of SO (resp. SI) (see Appendix A). 

4. Theoretical calculation of T for HD 49933 

We derive the mode excitation rates f for HD 49933. 
According to iGillon fc Magainl (|2006f l. HD 49933 has 
[Fe/H]=-0.37 ± 0.03 dex, while we only have two 3D sim- 
ulations with values of [Fe/H], respectively [Fe/H]=0 and 
[Fe/H]=-1. 

As seen in Sect. 13. 2[ differences in V between SO and 
SI are a direct consequence of the differences in the source 
functions Sr and Ss- It follows that in order to derive T 
for HD 49933, we only have to derive the expected val- 
ues for 5*^, and Ss- As seen in Sect. 13.21 differences in Sr 
(or in Ss) between SO and SI are related to the surface 
metal abundance through the surface densities that impact 
the convective velocities {u) . The determination of the HD 
49933 convective velocities allows us to determine its source 
function. To this end, we use the fact that the kinetic flux 
is almost unchanged between SI and SO (see Sect. 13. 2p to 
derive the profile of u(T), expected at the surface layers of 
HD 49933. This is performed by interpolating in Z between 
SO and SI, the surface density stratification representative 
of the surface layers of HD 49933. The whole procedure is 
described in Appendix A. 

In order to compute CP for HD 49933, we then need 
to know Z for this star. Since we do not know its sur- 
face helium abundance, we will assume by default th e 
solar value for Y- Y = 0.249 ± 0.003 (iBasiJ fl99l . 
iGillon fc Magainl ()2006[ )'s analysis shows that the chemi- 
cal mixture of HD 49933 does not significantly differ from 
that of the Sun. According to lAsplund et all (|2005| ). the 
new solar metal to hydrogen ratio is {Z/X)q = 0.0165 
Accordingly, since [Fc/H]=- 0.37 ± 0.03 dex, we derive 
Z = 5 3 10"^ ± 0.4 10'^ for HD 49933. Note that as- 
suming iGrevesse fc Noeld ()1993[) 's chemical mixture yields 
Z = 7.8 10-3 ± 0.5 10-3. 

The result of the calculation is shown in Fig. [21 
The maximum ? is 1.08 ± 0.05 lO^"^ J/s when 
lAsplund et al.l (|2005f )'s chemical composition is assumed 



(see Appendix A). This is about 30 times larger than in 
the S un and about 14 t imes larger than in a Cen A. 
When [G revesse fc Noelsl (|l993l )'s chemical mixture is as- 
sumed, the maximum in f is in that case equal to 
1.27 ± 0.05 ipl'^J/s, that is about 30 % larger than with 
lAsplund et al.l (|2005f )'s solar chemical mixture. 

We note that the uncertainties in the knowledge of 
[Fe/H] set uncertainties on !P which are the on order of 
10 % in the frequency domain of interest. 

5. Conclusion 

Wc have built two 3D hydrodynamical simulations repre- 
sentative in effective temperature (Toff ) and gravity (5) 
of the surface layers of an F type star on the main se- 
quence. One model has a solar iron-to-hydrogen abundance 
([Fe/H]=0) and the other has [Fe/H]=-1. Both models have 
the same Tcff and g. For each 3D simulation, wc have com- 
puted an associated "patched" ID full model. Finally, we 
have computed the mode excitation rates T associated with 
the two "patched" ID models. 

Mode excitation rates associated with the metal poor 
3D simulation are found to be about three times smaller 
than those associated with the 3D simulation which has a 
solar surface metal abundance. This is explained by the fol- 
lowing connections: the lower the metallicity, the lower the 
opacity. At fixed effective temperature and surface gravity, 
the lower the opacity, the denser the medium at a given op- 
tical depth. The higher the density, the smaller are the con- 
vective velocities to transport the same amount of energy 
by convection. Finally, smaller convective velocities result 
in a less efficient driving. On the other hand, a surface metal 
abundance higher than the solar metal abundance will re- 
sult in a lower surface density, which in turn will result in a 
higher convective velocity and then in a more efficient driv- 
ing. Our result can then be qualitatively generalised for any 
surface metal abundance. 

By taking into account the observed surface metal abun- 
dance of the star HD 49933 (i.e. [Fe/H] =-0.37), we have de- 
rived, using two 3D simulations and the interpolation pro- 
cedure developed here, the rates at which acoustic modes 
are expected to be excited by turbulent convection in the 
case of HD 49933. These excitation rates CP are found to be 
about two times smaller than for a model built assuming a 
solar metal abundance. These theoretical mode excitation 
rates will be used in Paper II to derive the expected mode 
amplitudes from measured mode linewidths. We will then 
be able to compare these amplitudes with those derived for 
HD 49933 from different seismic data. This will constitute 
an indirect test of our procedure which permits us to in- 
terpolate for any value of Z the mode excitation rates CP 
between two 3D simulations with different Z but the same 
Toff and logg. We must stress that a more direct valida- 
tion of this interpolation procedure will be to compute a 
third 3D model with the surface metal abundance of the 
star HD 49933 and to compare finally the mode excitation 
rates obtained here with the interpolation procedure with 
that obtained with this third 3D model. This represents a 
long term work since several months (about three to four 
months) are required for the calculation of this additional 
3D model, which is in progress. 
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Appendix A: Theoretical calculation of the mode 
excitation rates for HD 49933 

The mode excitation rate T is inversely proportional to the 
mode mass M (see Eqs. © and This is why we 
can derive M and MT separately in order to derive V for 
HD 49933. 

A.l. Derivation ofMT 

As pointed out in Sect. 13. 2i the kinetic flux Fkin ~ pv? 
is almost unchanged between SI and SO because both 3D 
models have the same Teff. This has also to be the case 
for HD 49933 (same T^s and same log g than SO and SI). 
Therefore, the calculation of M T for HD 49933 relies only 
on the evaluation of the values reached - at a fixed mode 
frequency - by the source functions §/{ and §s. 

As seen in Sect. 13. 2[ ujq = k^u controls the width of Xk 
in a way that the source functions §_r(wosc) and §s(wosc) 
can be seen as functions of the dimensionless ratio Wosc/'j-'o- 
The variation of E with k as well as the variation of Xk 
with lo/lji^) are shown to be similar in the two 3D simula- 
tions. Furthermore, SO and SI have approximately the same 
characteristic length A and hence approximately the same 
ko = 27r/A. Therefore, the source function (rcsp. Ss) 
associated with SO only differs from that of SI by the char- 
acteristic velocity u. This must then also be the case for 
HD 49933. Further, in order to evaluate the source func- 
tions in the case of HD 49933, we only need to know the 
factor 7 by which u is modified in HD 49933 with respect 
to SI or SO. According to Eq. ^ (rcsp. Eq. multipling 
u by 7 is equivalent to replace Sfl(ti'osc) (reps. §5(a;osc)) by 
7§flXwosc/7) (resp. 7 S,s(cJosc/7))- 

Since the kinetic fiux i^kin in HD 49933 must be the same 
for SO or SI, the characteristic velocity u can be derived 
for HD 49933 according to u*(T) = ui 7, with 7*(T) = 
{pi/p*y^^ where pi{T) is the mean density stratification of 
SI, ui{T) the characteritic velocity of SI and p* the mean 
density of HD 49933. Once 7* and then u^, are derived for 
HD 49933, we then compute the source functions associated 
with HD 49933. Finally, we compute M T by keeping Fkin 
constant. We now turn to the derivation of the factor 7*. 

A. 2. Derivation ofj^ 

To derive 7* at a given T, we need to know how the mean 
density p varies with the metal abundance Z. In order to 
this we consider five "standard" ID models with five differ- 
ent values of Z. These ID models are built using the same 
physics as described in Sect. 12.31 Two of these models have 
the same abundance as SO and SI. All of the ID models 
have approximately the same gravity (logg — 4.25) and 
the same effective temperature (Toff ~ 6730 K). 

The set of ID models shows that - at any given temper- 
ature within the excitation region - p varies with Z rather 
linearly. In order to derive p for HD 49933, we apply - at 
fixed T and between SO and SI - a linear interpolation of 
p{T) with respect to Z. 

A. 3. Derivation ofM 

As shown in Sect. 13.21 above in the frequency domain where 
modes are detected in HD 49933, M does not change sig- 
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nificantly between SO and SI. This suggests that the mode 
masses associated with a patched ID model with the metal 
abundance expected for HD 49933 would be very similar 
to those associated with SO or SI. Consequently we will 
assume for the case of HD 49933 the same mode masses 
as those associated with SI, since this 3D model has a Z 
abundance closer to that of HD 49933. 

A.4. Derivation of 7 

Before deriving T for HD 49933, we check that, from SO and 
the knowledge of ^Pq, we can approximately reproduce ^i, 
the mode excitation rates, associated with SI following the 

1/3 

procedure described above. Let 71 = (po/pi) . As seen in 
Fig. [3l when we multiply uq by 71 (T) we matche ui. Then, 
using 71 (T) and following the procedure described above, 
we derive Toi, the mode excitation rates associated with 
SI but derived from SO. The result is shown in Fig. lA.ll 
Tqi matches Ti rather well. However, there are differences 
remaining in particular in the frequency domain v=l.2- 
1.5 mHz. Nevertheless, the differences between CPqi and T 
are in any case not significant compared to the accuracy at 
which the mode amplitudes are measured with the CoRoT 
data (see Paper II). This validates the procedure, at least 
at the level of the current seismic precisions. 
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Fig. A.l. Mode excitation rates CP as a function of the 
mode frequency v. The thin dot-dashed line corresponds 
to Vqi, the mode excitation rates derived for SI from SO 
(see Appendix A4). The other lines have the same meaning 
as in Fig. O 



Since the metal abundance Z of HD 49933 is closer to 
that of SI than that of SO, we derive the mode excitation 
rates CP associated with HD 49933 from SI following the 
procedure detailled above. The result is shown in Fig. [2l 
As expected, the mode excitation rates 7 associated with 
HD 49933 lie between those of SO and SI, while remaining 
closer to SI than to SO. Note that the differences between Ti 
and the excitation rates derived for HD 49933 (CP) are of the 
same order as the differences seen locally between CPi and 
CpQi- These differences remain small compared to the cur- 
rent seismic precisions. On the other hand the differences 
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between 7 and Tq are significant and have an important 
impact on the mode amplitudes (see Paper II). 



